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In this contribution, we present some predictions for the production of D and B mesons in ultra- 
relativistic heavy ion collisions at RHIC and LHC energies and confront them with experimental 
results obtained so far by the STAR, PHENIX, ALICE and CMS collaborations. We next dis- 
cuss some preliminary results obtained with an improved description of the medium based on 
EPOS initial conditions, and its possible implications on the nuclear modification factor and on 
the elliptic flow of heavy quarks. 
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1. (Short) Introduction 

Since the discovery that heavy quarks (HQ) manifest a large coupling with the hot medium 
created in heavy ions ultrarelativistic collisions, many models and theoretical schemes (for a review 
see e.g. [JXJ, ||]) have been proposed in order to achieve a quantitative understanding of experimental 
data [g]. While a consensus has not been achieved yet on the most appropriate way to treat the 
coupling of heavy quarks with this medium, new results from LHC experiments provide more 
constraints on the various approaches and enrich the debate. In this contribution, we recall some 
of our results for RHIC experiments and provide the ensuing predictions for LHC energies, both 
obtained by taking the fluid dynamic approach of Kolb and Heinz (KH) to model the background 
medium. We then present some preliminary study of the consequences of taking an improved 
description of the medium by choosing the EPOS initial conditions. 

2. Results at RHIC 

In fig. [T[ we recall predictions from our "Monte Carlo at Heavy Quarks" generator (MC@. V HQ) 
for the nuclear modification factor of non-photonic single-electrons (NPSE) stemming from D and 
B decay in Au-Au collisions at yfs = 200 GeV / c using either collisional energy loss or collisional 
+ radiative energy loss for HQ as described in [Eh. Both models are able to cope with the NPSE 
Raa data measured by RHIC experiments, at the price of a global rescaling of the interaction rate 
(K = 2 and 0.65 respectively). As explained in [g], this model "fragility" has its origin in the rather 
stiff initial p t -distribution of heavy quarks, which makes a sequence of many small losses more 
probable than a single process involving large energy loss. 
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Figure 1: Raa of non photonic single electrons for central Au-Au collisions at RHIC energy and mid- 
rapidity; see Q| for details. 

One thus recovers Langevin's picture, so that the full dynamical evolution is encompassed 
in a single transport coefficient, e.g. the drag or the diffusion coefficient. Keeping K fixed and 
considering other centralities or the elliptic flow observable, one obtains a very good agreement as 
well, which seems to confirm this interpretation. In [Q], we have presented the drag coefficient A for 
the various rescaled models that allow to match RHIC NPSE data and have extracted values for the 
spacial diffusion coefficient D s of 2nTD s « 1.9 ±0.5 at T = 200 MeV and of 2nTD s « 2.6 ± 0.65 at 
T = 300 MeV, for both c and b quarks, in quite good agreement with recent lattice calculations [||]. 

In fig. |2], we show our predictions separately for the D and the B mesons for the same system as 
in fig. |I[ Accordingly, one sees that it might be difficult to distinguish between both types of energy 
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loss mechanisms, unless future upgrades of RHIC experiments permit to access heavy mesons with 
p t > 10 GeV/c. We have checked that the peak observed at p t 1.5 GeV/c in Raa(D) is essentially 
due to the radial flow of light-quarks that recombine with heavy ones through coalescence in order 
to form heavy mesons and not due to the radial flow of heavy quarks themselves (the Raa (c) is 
peaked at p t = GeV/c), which explains why such peak is not observed in Raa{B). 
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Figure 2: Prediction for the Raa of heavy mesons according to the same models as those used in fig. |jj 



3. Predictions at LHC 

We have made predictions [Q, [7|, |8|] for the production of heavy mesons in Pb+Pb collisions at 
^/s = 2.76 TeV using the same elementary cross sections as for the RHIC case and adapting the 
initial entropy density so in the fluid dynamics to obtain 



dy 



1600 at the chemical freeze out as 



well as the initial /^-distributions (taken from FONLL 1.3.2). Fig. [3j shows our latest calculations 
for the Raa of heavy mesons at mid-rapidity (recent predictions for the V2 and the centrality depen- 
dence can be found in [|, §]). The "El. + LPM" curve on the left panel corresponds to our model 
where coherence effects are taken into account for the radiative energy loss by considering a single 
effective scattering center; for B mesons, those effects have little consequences. 
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Figure 3: Same as fig. ||for y/s = 2.76 TeV; see text for details. 



The quenching of B mesons has not been measured directly yet but can be bona fide assimilated 
to the quenching of non-prompt J/\j/ measured by CMS with Raa ~ 0.35 for 6 GeV/c < 
p 3 J^ < 30 GeV/c, in good agreement with our predictions. For D mesons, our results show an 
excess of quenching as compared to the recent ALICE data [10] - especially when one includes a 
significant proportion of radiative energy loss - which is under present study. Interesting though is 
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the possibility to discriminate between collisional and radiative energy loss provided one acquires 
data in the p t range of 10-30 GeV/c. 



4. Influence of an improved medium description 



In the following we will discuss the influence of an improved description of the fluid dynamic 



background which we choose from EPOS initial conditions [11]. Including a final hadronic cas- 
cade this model is able to simultaneously describe a variety of soft observables such as particle 
yields, spectra, flow coefficients and dihadron correlations at RHIC and LHC energies. Having the 
soft sector under control gives us confidence that we can reliably investigate the medium-induced 
energy loss of heavy quarks. The features which distinguish the fluid dynamic evolution of the 
EPOS initial conditions from the previously used KH fluid dynamic background are explained in 
the following. The initial conditions are obtained from a flux tube approach to multiple scattering. 
Each elementary scattering process is described by a partem ladder, whose final state is a longitu- 
dinal color field. The dynamics of this flux tube is described by a relativistic string. In elementary 
collisions the string breaking by qq production leads to hadron formation from the individual string 
segments. In nucleus-nucleus collisions the density of flux tubes is large and string segments, 
which are slow and/or far from the surface, are transformed to fluid dynamically evolving matter. 
From this procedure one obtains the initial profiles for all fluid dynamic fields including the initial 
transverse velocities. These latter quantities turn out to be very small and can be adjusted as a 
parameter. This random flux tube distribution allows for an event-by-event treatment of the fluid 
dynamic evolution accounting for the fluctuating spatial structure of single events. In principle we 
could study the heavy quark propagation in a realistic event-by-event simulation. For the gain of 
statistics, however, we evolve several hundred thousand MC@sHQ runs per fluid dynamic event. 

The fluid dynamic evolution is governed decisively by the underlying equation of state (EOS). 
While the KH fluid dynamic medium is obtained from an unrealistic EOS with a strong first order 
phase transition based on ideal hadron and quark-gluon gases, subsequent to the EPOS initial con- 
dition the full 3 + Id fluid dynamic simulation is performed including a parametrization of the EOS 



from lattice QCD Q12|]. It exhibits a crossover transition between partonic and hadronic degrees of 
freedom in a range of temperatures between T = 145 — 165 MeV fill 
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Figure 4: Temperature evolution and radial velocity profiles for a central RHIC event comparing the EPOS 
to the KH medium. 
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The key ingredients for the propagation of heavy quarks are the temperature fields and the 
fluid velocities. These two quantities from one consistent fluid dynamic evolution are shown in fig. 
Q The time evolution of the temperature in the center of the collision area is shown in the left plot 
of fig. [| for the fluid dynamic evolution of EPOS compared to that of KH. The initial time for the 
fluid dynamic expansion is To = 0.35 fm in the case of EPOS and To = 0.6 fm for KH. A clear 
difference is seen at later times, when the KH EOS exhibits a first order phase transition leading 
to a constant temperature of 165 MeV, while the temperature in the crossover case keeps dropping. 
Although the temperature evolution above the phase transition is similar for both media, due to the 
different EOS the energy density, pressure and the resulting velocities are different. The right plot 
of fig. |] shows the radial velocity profile at different times comparing the EPOS to the KH medium. 
The radial velocity in EPOS is finite initially and then develops quickly toward its final value up to 
around 0.9c at the edges of the fluid. There is no initial radial velocity in the KH medium and the 
maximal values reached during the evolution are below those from EPOS. 

We initialize the heavy quarks spatially at the nucleon-nucleon collision points given by the 
elementary scatterings in EPOS in the transverse plan at zero space-time rapidity. Figure || demon- 
strates the influence of the initial radial velocity on the heavy quark observables: the nuclear mod- 
ification factor Raa (left plot) and the elliptic flow V2 (right plot). These simulations are done for 
Au+Au collisions at a beam energy of ^fs = 200 GeV and a centrality — 10% {Raa) or a fixed 
impact parameter b = 7 fm (V2), with the "El. + LPM" cocktail and K = 1. The Raa corresponds 
to the spectra of c-quarks at the transition temperature T = 155 MeV, where we stop the evolution. 
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Figure 5: The Raa an d the elliptic flow for charm quarks for different initial radial velocities at RHIC. 



In the Raa of c-quarks one can observe how the p t spectra are shifted toward higher p t for 

(2) 

increasing initial radial velocities, where ini corresponds to the standard value utilized in EPOS 
and for which the profiles were shown in figure 0. At the largest initial radial velocity investigated, 
v rad ini tne s P ectra at l° w Pt are depleted and one finds a peak in Raa at around p t — l GeV. For the 
investigated range of initial radial velocities the Raa is influenced in the low and intermediate p t 
range up to p t ~ 7 GeV. The quenching for high-/? f particles is obviously the same irrespectively 
of the initial radial velocity. A larger initial radial velocity accelerates the radial expansion of the 
fluid and counteracts the building up of the elliptic flow at small-intermediate p, [14]. This effect 
becomes more evident for particles with larger mass. Within our approach we can clearly observe 
this issue in the V2 of c-quarks, see right plot of figure |5[ For the largest initial radial velocity we 
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find the smallest V2 up to p t ~ 2.5 GeV, where V2 is largest for zero initial radial velocity. 

Finally, we compare the c-quark /?aa (10% most central) and V2 (fixed b = 7 fm) for the EPOS 
and the KH medium in the left and right plot of figure || respectively. For this purpose we stop the 
evolution at T = 168 MeV, which is above the first order phase transition in KH, and take K = 0.7. 
Due to the larger radial velocity in the EPOS medium the spectrum of c-quarks is shifted toward 
higher p, compared to KH. Therefore, the c-quark R^a atp t <l GeV is much lower in the EPOS 
medium than in the KH medium. If one switches off the radial velocity for the quark propagation 
by hand both media yield the same 7?aa at small p t . A slight difference remains stemming from the 
small differences in the temperature profiles, the fluctuating initial conditions and the different Xq. 
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Figure 6: The R^ and the elliptic flow V2 of charm quarks for the EPOS and the KH medium with and 
without radial velocity, central RHIC events. 



If the radial velocities are switched off completely the only contribution to the elliptic flow is 
p t broadening over the different path lengths in x and y directions. Therefore, we find a reduced 
but finite elliptic flow V2 which is of the same order in the EPOS and the KH medium. Since the 
EOS are different, the flow pattern in EPOS and KH yield significantly different elliptic flow if the 
full radial velocities are considered. 
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